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AN ANALYSIS OF THE EFFECT OF DESIGN AND MEASUREMENT ERRORS 
ON PUMP PERFORMANCE PARAMETERS 
by Douglas A. Anderson 
Lewis Research Center 

SUMMARY 

A method is presented for calculating the effects of design and measurement errors 
on pump performance parameters. Error equations and charts are presented which 
relate the amount of error in a given performance parameter to the amount of error in 
a given design or measured variable. The design error equations were developed for a 
blade element design procedure. The measurement error equations apply to pump tests 
in which radial surveys of pressure and angle measurements are made. 

Some specific applications of the error analysis equations are discussed. The com- 
puted error values provide a criterion of how large a change in a performance parame- 
ter must be to be considered significant. This is useful when comparing design per- 
formance to measured performance or when comparing the measured performance of 
different pumps. The error equations may also be used to determine the sensitivity of 
a given design to design errors and manufacturing tolerances, or to assess the relative 
importance of several measurement errors in a particular test configuration. The ap- 
plications often require that the effect of two or more errors on a given performance 
parameter be considered. This is done by using the concept of the uncertainty interval, 
which is explained in this report. Finally, an example is given which demonstrates the 
application of the measurement error analysis to some typical test data. 

Some general trends of the sensitivity of the performance parameters to design and 
measurement errors were deduced. The sensitivity to design errors generally in- 
creased with increasing stagger angle or decreasing loading. For measurement errors, 
no general trends with stagger angle or loading were indicated. 

The error equations were developed for axial flow rotors, but are not limited to 
this type of blade row. Modifications to the equations are presented which allow them 
to be applied to stators. Another modification is mentioned which allows them to be ap- 
plied to mixed-flow impellers. 



INTRODUCTION 


The importance of measurement errors in a pump test program is generally rec- 
ognized. In fact, it would be hard to evaluate any set of test data without developing 
some intuitive feeling about the effect of errors. However, mathematical techniques are 
available which allow a methodical analysis of the effect of measurement errors on the 
resulting flow and performance parameters. Using these techniques, the amount of 
error to be expected in the performance parameters can be calculated from estimates of 
the errors introduced by the different measurements. An error analysis can also be of 
use in planning instrumentation by showing which measurement is likely to introduce the 
largest errors. 

An analysis of the effect of errors should also be useful in a pump design system. 
Such an analysis would show the sensitivity of a given design to design errors and manu- 
facturing tolerances. It would allow estimates of design precision, the limits around the 
design point within which the pump performance may be expected to lie with reasonable 
certainty. Also, an error analysis would show which step in a design procedure is the 
most critical source of error. These examples show just a few of the applications of 
an error analysis to pump design procedures and test results. 

An error analysis was developed in reference 1 and applied to axial flow compres- 
sors. Equations were derived which relate the amount of error in a given compressor 
performance parameter to the amount of error in a selected design variable or measured 
variable. The methods used herein to derive the error equations for axial flow pumps 
are the same as used in reference 1, but almost all of the derived error equations are 
different. They are different because different parameters are used to describe pump 
performance, due to the incompressibility of the pump flow, and the effect of cavitation. 
Further differences between the pump and compressor error analyses arise because of 
different ranges of interest for some parameters (e. g. , blade stagger angle and flow 
coefficient) and slightly different measurement techniques. 

The purpose of this report is to present the methods which can be used to calculate 
the effects of both design errors and measurement errors on pump performance. Equa- 
tions are derived which relate the amount of error in a given pump performance param- 
eter to the amount of error assumed in a given design or measured variable. These 
equations apply to flow along streamlines or along blade elements. The equations are 
presented in the form of carpet plots. Whenever feasible, the pump error analysis is 
presented in a manner similar to that of reference 1. Some uses of the error analysis 
equations are suggested. The effects of changes in stagger angle and loading on the 
sensitivity of the performance parameters to design and measurement errors are dis- 
cussed. An example is given in which the measurement error equations are applied to 
some typical pump test data. 
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DERIVATION OF EQUATIONS 


General Error Analysis Equations 

Consider a function of n variables, f(xp x 2 , . . . , x n ). The change in f due to 
a small change in one of the variables x m is given approximately by 


Af m = 


af 


Ax 


3x. 


m 


m 


(i) 


(Symbols are defined in appendix A. ) If the error in each variable x^, Xg, . . . , x n is 
known, the resulting error in the function f can be computed from 


Af = -§L. Ax, + AXr, + . . . + Ax (2) 

0Xj 3xg 9x n 

where A denotes a known or calculated value of the respective error. This equation 
can be used for systematic errors, which are errors of known sign and approximately 
known magnitude. 

However, the errors encountered in the design and measurement analyses are more 
likely to be of the type called random errors or uncertainties. Reference 2 discusses 
the nature of these uncertainties and how they should be treated statistically. (Uncer- 
tainty, as defined in ref. 2, is the more precise term. ) The relative size of the error 
or uncertainty can be given in terms of an uncertainty interval as used in references 
1 and 2. If, for example, for a measured pressure of 50 psi, there is an estimated 
90 percent probability that the true pressure is between 49 and 51 psi, the uncertainty 
interval is 2 psi for this probability. If the uncertainty intervals of the other measured 
variables at the same probability are known, the uncertainty interval of the parameter 
f for this probability can be computed from 



In this equation, A refers to an uncertainty interval. (Except for a difference in nota- 
tion, this is the same as eq. (7) of ref. 2.) It is assumed that all of the errors AXj, 
AXg, etc. , have a normal distribution, but the equation is sufficiently accurate for any 
distribution likely to occur. It is also assumed that the variables x^, Xg, etc. , are in- 
dependent; that is, a change in one will have no effect on the values of the others. 
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It is often more convenient to d'teal with dimensionless terms Af/f, Ax^/x^, etc. 
These terms are easily expressed as percentage errors, and a simplified form of the 
error equations often results. Using dimensionless terms, equation (2) for systematic 
errors might become 


Af /i_9L\ A x 2 + 

f \t a xjx l \t dx 2 J 


( 4 ) 


This equation shows the error Axj/xj in dimensionless form, but not Ax 2 - The mod- 
ified error equations will generally have a mixture of both forms. The form Ax is 
better for angles, but the dimensionless form, Ax/x, is usually more convenient for all 
other variables. The corresponding equation for random errors (eq. (3) modified) is 



X 1 3f 
f 3Xj 



+ 



(5) 


Design Errors 

The design error analysis was developed for a blade element design procedure such 
as is discussed in reference 3. In a blade element design procedure, the hydrodynamic 
design (calculation of velocity diagrams and performance parameters) is conducted along 
selected blade elements. Blade sections to achieve the desired performance are then 
selected. In a design procedure, certain quantities, which shall be called the design 
variables, are chosen by the designer or determined from appropriate design correla- 
tions. The design variables selected for this analysis are 
Inlet and outlet radius, r^ and 
Inlet and outlet axial velocity, and V ^ 

Wheel speed at the tip, U t 
Inlet total head, H.^ 

Inlet flow angle, 

Outlet relative flow angle or deviation angle or 6) 

Loss coefficient, co 

From the design variables, the pump performance parameters are computed using 
the design equations (eqs. (Bl) to (B12), appendix B). The performance parameters 
included in the design error analysis are 



Inlet and outlet velocities, and Vg 

Inlet relative flow angle or incidence angle, (3^ or i 

Outlet flow angle, 

Ideal head-rise coefficient, if/^ 

Head-rise coefficient, i]/ 

Blade element efficiency, 77 
Diffusion factor, D 

Although AHj and AH are often used instead of t//. and \j/, they are not included 
among the performance parameters, since the percentage errors are the same for AH. 
and and for AH and that is, A(AH i )/AH i = AiJ/./ty/j and A(AH)/AH = A^//^/. 

The values of the design variables, together with their inherent inaccuracies, are 
the necessary inputs required to do an error analysis of a design procedure. The values 
of the design variables v zi , p v and Hjl are known for selected values of r ^ from the 
analysis of the flow through the previous blade row or, for the first stage rotor, from 
the analysis of the flow through the pump inlet. The errors in V z ^, j 3 j, and H^ are 
the result of the accumulation of errors in all upstream blade rows. To estimate the 
errors in the design variables 6, c o, rg, and V Z 2, which are obtained from correlations 
with other parameters (e. g. , the ZD against D correlation, radial equilibrium equa- 
tion, etc. ) requires consideration of both the errors in these other parameters and the 
accuracy of the correlations used. As blade elements should ideally lie along stream- 
lines, any difference between an assumed or calculated rg value and that of a stream- 
line through the corresponding inlet radius is regarded as an error in rg. 

The design variable Hj is not included in the design error formulas, since none 
of the performance parameters depend directly on H^. No equations are given for de- 
sign errors in r^, as r ^ is regarded as a selected, nonvarying input. Errors in U t 
are regarded as measurement errors rather than design errors and consequently are 
not included in the design error formulas. However, a designer may wish to include 
their effect in estimates of design precision (how closely the actual performance is ex- 
pected to meet the design specifications). 

For systematic design errors, the error in a particular performance parameter f 
is given by 

" = (hi _xYhi + (1 jl W + (hi Ji-Yhi + (1 W 

£ \ f wzihzi v wj V ‘ w z 2 )v z2 \tse'J 

+ (1 iVw + /i? ( 6 ] 

V V 9r 2/ r 2 
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This is simply equation (4) with Xj, x 2 , etc. , identified as appropriate design variables. 
For the case of random errors, the uncertainty interval in f is related to the uncer- 
tainty intervals in the design variables by 




+ 



(7) 


These equations were derived assuming that the design variables are independent. 
However, in most design systems some of the design variables are interrelated (w, D 
correlation, etc. ) so that an error in one may provide errors in some of the others. 

The analysis of such interactions is beyond the scope of this report. However, the ef- 
fects of these interactions should be small compared with the direct effects of design 
errors with one exception. The outlet axial velocity V z2 is strongly affected by the 
values of |3 2 and through the radial equilibrium equation. 

The treatment of the variable r 2 requires careful consideration. Like all of the 
other variables considered in the error analysis, r 2 will have a certain amount of error. 
But as the outlet design and performance parameters are usually specified or calculated 
at given values of r 2 , it is not convenient to have to deal with errors in r 2 . Errors in 
r 2 are especially troublesome if the calculations are to be extended downstream to sub- 
sequent blade rows. One remedy is to modify the error values assigned to the other 
parameters to include the effects of uncertainties in r 2> Then r 2 can be treated as if 
its values were fixed. 

Assume that a streamline, nominally at radius r 2 , it actually at radius r 2 + Ar 2> 
The error in f at r 2 + Ar 2 would be given by the appropriate equation (eqs. (2) to (7)) 
except that the Ar 2 term would be zero. Since V z2 , /3 2 , and u> are the only design 
variables affected by a change in r 2 , the value of f at r 2 is related to the value of f 
at r 2 + Ar 2 by 


f(r 2 ) = f(r 2 + Ar 2 ) - 


3f dV Z2 f 5f d/3 2 
a v z2 dr 2 d(3' 2 dr 2 


3co dr 2 y 


( 8 ) 


where dV z2 /dr 2 simply refers to the slope of the graph of V z2 plotted against r 2 , 
and similarly for |3 2 , c o, and f. From the definition of the total derivative, 
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df _ 9f t 9f dV Z2 ( 3f d<3 2 t at dco 
dr 2 3r 2 aV Z2 dr 2 9 ^2 dr 2 9 “ ^2 

Thus, the following equation is equivalent to equation (8), 

f(r„) = f(r 2 + Ar 2 ) + f— - — W 2 

2 2 2 W 2 dr 2 ; 

Consequently, to get the correct value of Af (if the value of r 2 is to be regarded as 
fixed), 9f/3r 2 in equations (2) to (7) should be replaced by 

9f dV Z2 + J)f_^2 | 9f dco \ 

9V Z2 dr 2 9 ^2 dr 2 9 ^ dr 2/ 

or the equivalent expression 

9f df 
9r 2 dr 2 

The resulting error equations are 

(1) For systematic errors and a fixed r 2 value (using dimensionless terms as in 
eq. (6)), 



( 9 ) 


( 10 ) 
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The replacement of 3f/3r 2 by [( 3f/3r 2 ) - (df/dr 2 )] should not be considered as a sub- 
stitution, and the two expressions are not equivalent. Equations (6) and (7) are for an 
r 2 which has an error Ar 2 while equations (11) and (12) are correct for a fixed r 2 . 


Measurement Errors 

The measurement error analysis applies to a pump test in which radial surveys of 
total and static pressure and of flow angle measurements are made, as was done in ref- 
erences 4 and 5. The measured variables considered in this analysis are 
Inlet and outlet velocity, Vj and V 2 
Inlet and outlet flow angle, /3^ and /3 2 
Inlet and outlet radius, r^ and r 2 
Net positive suction head, H_„ 

bv 

Head rise, AH 
Tip speed, U t 

Velocity is not a directly measured quantity, but is computed from the measured 
value of velocity head. However, it is more convenient to use velocity in the error equa- 
tions, and the percentage error in velocity is simply half the percentage error in veloc- 
ity head. When velocity head is measured directly, as done in references 4 and 5, the 
accuracy is much better than if the velocity head were obtained from separate measure- 
ments of total pressure and static pressure. 

From the measured variables, the pump performance parameters are calculated 
from appropriate equations (eqs. (B46) to (B56)). The performance parameters included 
in this analysis are 

Inlet and outlet flow coefficient, <p ^ and cp 2 
Inlet and outlet relative velocity, Vj and V 2 
Inlet relative flow angle or incidence angle, /3^ or i 
Outlet relative flow angle or deviation angle, 0 2 or 6 
Head-rise coefficient, x(/ 

Ideal head-rise coefficient, \f/^ 

Blade element efficiency, r] 

Loss coefficient, w 
Diffusion factor, D 
Cavitation number, k 

The variables V z ^, V z2 , AH, and AHj are not included in the performance param- 
eters, since their percentage errors are the same as for cpp <p 2 , i ]/, and respec- 
tively, except when caused by errors in U^.. An error in U^. does not affect the accur- 
acy of the values of V Z1’ V Z2’ and AH, but the error in AH^ is given by 
A(AH i )/AH i = -A 
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The measurement error equations (eqs. (B57) to (B93)) relate the amount of error 
in a given performance parameter to the amount of error in a given measured variable. 
They are derived from equations (B46) to (B56) in appendix B. 

There are a number of factors which cause measurements errors. One is the in- 
herent errors of the instrumentation and data recording systems. The error estimates 
given in many data reports consider only this source of error. However, in addition to 
this, the effects of pressure and flow fluctuations, cavitation, etc. , on the pressure and 
angle measurements must be considered. Also, the values of r^ and rg may be in 
error due to inaccuracies in probe positioning. 

Since blade element test data are generally presented as plots of the measured 
variables and performance parameters against r ^ or rg, it would be preferable not to 
have any uncertainties in these variables. This can be done if the errors in the other 
parameters are modified to account for the errors in r^ and rg, using the method de- 
veloped in the design error section of this report. The error in a calculated perform- 
ance parameter f due to an error in radial probe position is given by 


where r stands for either r^ or rg. The total derivative df/dr is simply the slope 
of the curve of f plotted against r. The resulting error equations are 

(1) For systematic measurement errors, where r^ and rg are regarded as fixed 
values 



(2) For random measurement errors, where r^ and rg are regarded as fixed 
values 
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In the measurement error equations presented for Af/Ar in appendix B, the form is 
used which is judged most convenient for calculations. However, since r^ and rg are 
usually changed together, it is not possible to determine the individual values of df/dr ^ 
and df/drg from a graphically determined value of df/dr. The value of df/dr maybe 
used for df/dr g if f is not a function of r^ or for df/dr ^ if f is not a function of 
rg. For radially constant inlet total pressure, which is a normal test condition for inlet 
stages, dHg/drg can be replaced by dAH/drg. 


Presentation of Equations 

The error equations and the equations from which they are derived are listed in ap- 
pendix B. The equations are presented in the form judged most convenient for use in 
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TABLE I. - CLASSIFICATION OF THE ERROR EQUATIONS AND CARPET PLOTS 
ACCORDING TO PERFORMANCE PARAMETER AND DESIGN OR 
MEASURED VARIABLE 


(a) Design errors 


Design per- 





Design variables 





formance 

parameters 

V Z1 


h 

V Z2 

j3g or 6 

U) 


r 2 


Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Eq. 

Fig. 

v i 

B13 

— - 

B14 

1 

— 

— 






V 2 

— 

— 

--- 


B15 

2(a) 

B16 

2(b) 

— 

B17 

2(c) 

*1 

B18 

3(a) 

B19 

3(b) 


— 

— 


--- 

— - 



— 

--- 

— 


B20 

3(a) 

B21 

3(b) 

--- 

B22 

3(a) 

+ i 

B23 

... 

B25 

4(a) 

B26 

4(b) 

B28 

4(c) 


B29 

4(d) 

* 

B30 

5(a) 

B31 

4(a), 5(b) 

B32 

4(b) 

B32 

4(c) 

B33 

B32 

4(d) 

V 

B34 

2(a) 

B35 

2(b), 4(a) 

B36 

4(b) 

B36 

4(c) 

B37 

B36 

4(d) 

D 

B38 

6(a) 

B3° 

6(b) 

B40 

6(c) 

B41 

6(d) 

— 

B42 j 

7(a), (b) 


(b) Measurement errors 


Measured per- 


Measured variables 


>rmance pa- 

















rameter 


v i 

0 

1 

V 

2 

13 

2 

r 

1 

r 

2 

U t 

AH 

H sv 


Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Fig. 

Eq. 

Eq. 


B57 


B58 

1 









B59 

















w 2 

— 


--- 


B57 

--- 

B58 

1 

--- 

— 

--- 

— 

B59 

— 

— 

--- 

/3i or i 

B60 

3(a) 

B61 

2(a) 


— 

--- 


B62 

3(a) 

— 

— 

B63 

3(a) 

--- 

--- 

^2 or 6 

— 


--- 

--- 

B60 

3(a) 

B61 

2(a) 

— 

--- 

B62 

3(a) 

B63 

3(a) 

— - 

— 

* 

— 


— 




--- 

— 

... 


— 

— 

... 

B64 


B65 


*i 

B66 j 


B67 

— 

B6 r 


B69 

... 

B70 ! 

i 

4(a) 

B71 

4(a) 

B72 

... , 

— 


77 

B73 


B73 

— 

B73 


B73 

... 

B73 

4(a) 

B73 

4(a) 

B74 

— 

B75 

--- 

to 

B76 

2(a), 3(b) 

B77 

8(b) 

B78 

8(a) 

B79 ' 

8(b) 

B80 

--- 

B81 

... 

B82 

2(c) 

B83 

--- 

D 

B84 

6(a) 

B85 

9(a) 

B86 

9(b) 

B87 

9(c) 

B88 

--- 

B89 

... 

B90 

2(c) 

— 


k 

B91 

2(a), 3(b) 

B92 

8(b) 

— 

— 


— 

B93 

2(c) 

— 

— 

— 

— 

— 

B94 
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error calculations. Usually the equations are simpler if dimensionless terms such as 
(AVg/Vg) and [A?y/(1 - 77 )] are used. Some of the equations become much simpler when 
j3j is zero; when such is the case, the simplified form is included in the list of equa- 
tions. 

Carpet plots of some of the equations are presented in figures 1 to 9 to allow quick 
calculations. A carpet plot is simply a convenient way of showing a function of two 
variables. Each curve shown gives the relation between the function and one of the 
variables for a constant value of the other variable. For a function of more than two 
variables, a lattice arrangement of carpet plots can be used. Figure 6 (a) is an example 
showing a function of four variables. The individual carpets are offset so that a french 
curve may be used for interpolation between plots. For many of the error equations, 
two plots are presented, one plot covering the range of conventional pump operation, the 
other covering the higher blade angles typical of inducer operation. Most of the error 
functions had considerably different values and/or slopes at high blade angles as com- 
pared with moderate blade angles, so a single plot would not be readable at both ex- 
tremes. 

Not all of the error equations are presented in plot form, since for simple linear 
relations this was deemed unnecessary. Some plots are used for more than one error 
equation but, to avoid confusion, only one equation is represented in the plot coordi- 
nates. Table I gives the equation numbers and figure numbers for the various combina- 
tions of performance parameters and design or measured variables. Some of the plotted 
values must be changed in sign or multiplied by a constant in order to correspond to a 
particular measurement equation; for this reason, one should always refer to the error 
equations first before using the plots. 

The error equations were derived assuming small errors, but should be sufficiently 
accurate for most applications. An error equation gives the sensitivity to an error (i. e. , 
the ratio of the change in a performance parameter to the change in a design or meas- 
ured variable) at one particular set of conditions (i. e. , for particular values of /3j, 02 > 

<p g, etc.). A large error may change these conditions enough to significantly change the 
sensitivity. In this case, using the initial sensitivity to compute an error is not a good 
approximation. This possibility can be checked by comparing the sensitivity at two con- 
ditions differing by the amount of the assumed error. 


APPLICATION OF THE ERROR ANALYSIS 
Design Errors 

The design error equations and carpet plots (figs. 1 to 9) relate the magnitude of 


12 



the errors in selected performance parameters to the amounts of errors in individual 
design variables. As an example, take the highly loaded rotor data discussed in appen- 
dix C, where at a radius ratio of 0. 947 the measured values include i//. = 0. 391, 

/3g = 41. 7°, and |3 2 = 49. 0°. From figure 4(c), a 1° error in /3 2 results in an error 
in and \p of -4. 54 percent of the value of Since, for no prewhirl \p^ = a 

1° error in j 3 2 gives an error of -0.018 in \p^ and ip. Errors in 0 2 can be caused 
either by design errors in deviation angle or by inaccuracies in the fabrication of the 
blade outlet angle. The errors shown by the other design error equations can be deter- 
mined in a similar manner. 

The error equations can be used to determine the sensitivity of a given hydrody- 
namic design to design errors and manufacturing tolerances. The ratio of a 0. 018 error 
in ij/ to a 1° error in |3 2 may be regarded as a measure of the sensitivity of ip to 
errors in 0 2 . The sensitivity to errors may be important in comparing two competing 
designs. If it can be assumed that the size of the errors will not differ much in the 
designs being compared, the error equations will show which design is likely to meet its 
specifications the closest. It is not necessary to know the actual size of the errors in 
making this comparison. 

In an iterative design procedure, the output values of one iteration are generally 
used to determine the input to the next iteration. A typical design procedure might have 
each iteration begin with new values of 6, and r 2 . Using the design error equations 
or plots, it can be determined whether the difference between the old and new values of 
these inputs are great enough to significantly affect the other parameters. Another ap- 
proach is to terminate calculations when the differences between the new and old values 
of to and 6 are less than the anticipated design errors. Comparing the new and old 
values of r 2 is usually unnecessary. Values of r 2 are generally determined from a 
radial equilibrium calculation in such a way that specifying to and 6 will determine 
r 2 for a given geometry. To summarize, the error equations can be used to set defi- 
nite criteria on when to terminate the iterations. 

In comparing design performance to measured performance, an error analysis may 
help determine the significance of any discrepancy. The question is whether an observed 
discrepancy is within the limits to be expected for a given error or combination of 
errors. The combined effect of errors in two or more design variables may be calcu- 
lated using equation (2) or (3). A similar calculation may be made for measurement 
errors, or for a combination of design errors and measurement errors. 

Estimates of uncertainty intervals in the performance parameters are obtained if 
the effects of random errors in all of the design variables are combined using equa- 
tion (7). The uncertainty interval shows the range within which the parameter may be 
assumed to lie with reasonable probability. For example, a head-rise coefficient of 
0. 37 with an uncertainty interval of 0. 02 for 90 percent probability means that there is 
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a 90 percent chance that the value of i// is between 0. 36 and 0. 38. The uncertainty in- 
terval is the most reliable criterion of the accuracy of a set of data and is to be pre- 
ferred in making such comparisons as design performance with measured performance. 
However, the uncertainty interval is dependent on estimates of errors in all of the design 
variables. If one or more of the error estimates are questionable, it may be better to 
restrict the comparisons to the design parameters with reasonably accurate error 
values. 


Measurement Errors 

The measurement error equations and charts give the amount of error in selected 
performance parameters as functions of the amounts of measurement errors. One use 
of the measurement error analysis is to help in planning instrumentation. The error 
equations can be used to determine which measurement errors cause the greatest errors 
in the various performance parameters. Significant effort would be justified in order to 
improve the accuracy of a critical measurement, while only minimal effort should be 
applied to a measurement which introduces relatively small errors. 

The effects of random errors due to different measurements may be combined using 
equation (3). If the errors due to all pertinent measurements are included, the result 
is an uncertainty interval in the measured performance parameter. Since it includes 
the effects of all known random errors, the uncertainty interval is the best criterion of 
normal data variability. 

The measurement error calculation is useful in evaluating data trends. It gives an 
indication of the amount of data scatter which may reasonably be expected. If the scatter 
is small compared with the predicted amount of error, this indicates that the measure- 
ments are more consistent than anticipated. The error analysis may also show whether 
or not a change in the data is significant. For example, suppose that after a rotor is 
modified, the test data show a small increase in the values of blade element efficiency 
over the previous test data. It cannot be concluded that the modification improved the 
performance unless the differences are greater than can be explained by normal meas- 
urement errors. The amount of error to be expected is something to consider when 
evaluating any trend in a set of data. 

The measurement error analysis is applied to two particular sets of test data in 
appendix C. The data chosen for this example (refs. 4 and 5) were for a highly loaded 
conventional rotor and for a flat-plate inducer. The results are presented in terms of 
uncertainty intervals. Many of the ideas discussed in a general way in this report are 
brought out in this example, among them the trends of sensitivity with blade angle and 
the use of the error analysis in choosing instrumentation. The example is intended to 
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show some of the results which can be obtained from the error analysis and some of the 
problems involved in its use and to give an idea of what degree of accuracy to expect in 
current pump test data. 


TRENDS SHOWN BY ERROR ANALYSIS 

Using the carpet plots, some generalizations can be made on how the sensitivity to 
design and measurement errors will vary with different designs. The sensitivity of the 
performance parameters to design and measurement errors can be read directly from 
the curves of figures 1 to 9. The following paragraphs summarize the trends of sensi- 
tivity with changes in stagger angle and blade loading. 

The effects of a change in stagger angle on the sensitivities to design and measure- 
ment errors are listed in table II. The table shows how the sensitivity of a particular 
performance parameter (shown in the column on the left) to errors in a given variable 
changes when stagger angle is increased. The intention of this table is to allow general 
trends with stagger angle to be visualized. So if a particular trend differs depending on 
the values of the flow angles, no attempt is made to show the trend under all possible de- 
sign conditions. Rather, the trend under the most typical design condition is listed in the 
table, along with the restriction under which this trend applies. 

For design errors, (table 11(a)) the sensitivity trends show a strikingly regular pat- 
tern. The sensitivity of all of the pertinent performance parameters to errors in the 
outlet velocity diagram variables V Z 2 > anc * r 2 i ncrease as stagger angle in in- 
creased. (The sensitivity of D to errors in ^ is a possible exception. ) In constrast, 
the sensitivity to errors in the inlet velocity diagram variables and 0^ either de- 

crease or do not change as stagger angle is increased. For typical pump designs, the 
sensitivity to errors in V Z 2 > @ 2 > and ^ is considerably greater than the sensitivity to 
errors in V Z i and 0^. Consequently, the total amount of error in the performance pa- 
rameters follows the trend of the former variables. 

The effects of a change in blade loading on the sensitivity to design errors were also 
determined. (These results are not included in table II. ) The sensitivity of all of the 
performance parameters to errors in V Z 2 > £ 2 ’ an< * r 2 decreases as loading is in- 
creased. As loading has no direct effect on the inlet velocity diagrams, the sensitivities 
to errors in and 0^ are not affected. The general conclusion is that the sensi- 

tivity to design errors in all of the performance parameters except Vj and 0^ tend to 
increase with increasing stagger angle or decreasing loading. For an inducer, the com- 
bination of high stagger angle and light loading make the design very susceptible to de- 
sign errors. 

The sensitivity of the measured performance parameters to measurement errors 
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TABLE II. - TRENDS OF SENSITIVITY TO ERRORS WITH INCREASING STAGGER ANGLE 


(a) Design errors 


Perform- 



Design variable 



ance pa- 
rameter 

V Z1 

h 

V Z2 

4 

w 

r 2 


No change a 

NCT b 





v i 

<>i 

■ 




Decreases (high 
blade angle) c 

Decreases 

Increases (high 
blade angle) 

Increases 


Increases 

V 2 






^2 

Increases (ex- 
cept high /3g) 

Increases 
(02 ^ 02^ 

Increases (ex- 
cept high /3 2 ) 

No change 

i 

NCT 


* i 

Increases 
(02 > /3 2 ) 

Increases (high 
blade angle) 

Increases 

(^2 > ^ 2 ) 


* 

Decreases slightly 
(small Pj) 

Decreases slightly 
(small (3^) 

Increases 
(^2 ^ @2* 

Increases (high 
blade angle) 

Decreases 

slightly 

Increases 
(02 > 02^ 

V 

Decreases (small 

W 

Decreases (small 



NCT 

Usually decreases 
(small /}]) 

Increases 

<4 > **2> 

Increases 

Increases (high 
blade angle) 

Increases 

No change 

Increases 

<4 > 4 > 

NCT b 

D 



(b) Measurement errors 


Perform- 

Measured variables 

ance pa- 
rameter 

v i 


V 2 

02 

r i 

r 2 

u t 

AH 

*2 



No change 3- 

Increases 



No change 





NCT b 


4 

Decreases (high 
blade angle) c 1 

Increases (high 
blade angle) 

Decreases (high 
blade angle) 




No change 

* 




No change 

No change 

Decreases 

No change 

Decreases 

NCT 

NCT 


No change 



V 

No change 

Decreases 

No change 

Decreases 

NCT 

NCT 

No change 

No change 

w 

D 

Usually de- 
creases 

Usually de- 
1 creases 

Decreases (high 
blade angle) 

Usually de- 
creases 

, Usually in- 
creases 

Increases (high 
blade angle) 

Decreases (high 
blade angle) 

Usually de- 
creases 

NCT 

NCT 

[ 

NCT 

NCT 

NCT 

Increases 

Slight increase 



a Stagger angle does not change the sensitivity. 
b No consistent trend of sensitivity with stagger angle. 

c Parentheticals show conditions under which the trend applies. The restriction to high blade angles is never more severe than 
> 45° or /3 2 + /3g > 90°. 
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is also affected by changes in stagger angle and blade loading. However, unlike the de- 
sign error case, the sensitivity trends of the measurement errors showed no regular 
pattern. The individual trends with stagger angle are listed in table 11(b). For increases 
in blade loading, the sensitivity to errors either increases or does not change for all of 
the measurement error equations. The overall level of measurement errors (due to all 
sources of error), as compared with design errors, will in general be less affected by 
changes in stagger angle or blade loading. 

Some of the sensitivity trends are not immediately obvious from the carpet plots. 

For example, the sensitivity of to design errors in 0jj (fig. 4(c)) increases with 
incr easing 0lj, but decreases with increasing jSg* The trend with stagger angle will 
depend on which effect is stronger. In preparing table II, sample calculations were 
made for each case in which the trends in ^ an< * P 2 were opposed. It was found that 
whenever 0g was greater than 0g, the trend with 0g predominated. For most blade 
elements in a conventional pump design, 0^ is greater than 0g; consequently, the trend 
with jSg determines the trend of sensitivity with stagger angle for most blade elements. 

Application of Error Equations to Stators 

The design equations and measurement equations for stators differ from the related 

rotor equations in that, since U„ = 0, there is no distinction between relative and abso- 

s 

lute velocities, and \J/^ is zero. The stator error equations, which are derived from 
the stator design and measurement equations, also differ from the related rotor error 
equations. The two sets of error equations are, however, enough alike so that a few 
substitutions will produce the stator error equations from the rotor error equations. 
Therefore, even through the error equations presented herein were developed for rotors, 
they can, with suitable modifications, be applied to stators. 

The error equations for stators can be obtained from the listed error equations by 
making the following substitutions: 

(1) Whenever a relative angle appears, substitute the negative of the absolute angle; 
that is, substitute -0^ s for j8j, -02s ^ or @2’ anc * -A ^2s ^ or ^2' 

(2) For \p substitute 0. 

(3) For all other parameters, substitute the corresponding stator parameter. 

(4) In the error equations for D, change the sign of the second term (the term con- 

V 2 AV. 

tinuing a). This is required because, in the definition of D = 1 - — - + , the turn- 

Vj 2crV 1 

ing (A Y q term) is in the opposite direction for a stator than for a rotor. 

(5) These substitutions give the correct stator error equation with one exception, 
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equation (B42). The sign on the whole equation (not just the second term) should be 

changed and (r2V 0 2 “ r i V 0i)/^ r z2 r t substituted for 

After these substitutions are made, the resulting equation may simplify. For ex- 
ample, equation (B15) becomes 



The error equations for the performance parameters t], and become 
either meaningless or trivial when applied to stators and would not be included in a 
stator error analysis. For stators, rj is not defined, if/ i is identically zero, and j3 
and p '2 are the same as /3^ and /3 2 which are input variables in the stator error equa- 
tions. However, the error equations in which ^ appears as a design variable should 
not be dropped from the analysis. Using the substitutions, these equations are trans- 
formed into stator error equations in j 3 2 which are necessary for a complete error 
analysis. 

The stator counterpart of head-rise coefficient \p s can be defined using the value 
of from the preceeding rotor. The error equations derived for i£/ g will then give 
the effect of stator errors on the stage head-rise coefficient. 

For stators, the design and measurement equations for the inlet parameters V^ g 
and <p^ g are identical in form to those for the outlet parameters V 2g and cp 2g . As 
would be expected, the derived error equations show the same similarity. The inlet 
error equations are easier to derive and can then be applied to both inlet and outlet 
parameters. 


CONCLUDING REMARKS 

The basic purpose of this report is to provide a method for calculating the effects of 
design errors and measurement errors on pump performance. Error equations were 
developed which allow the amount of error in a given performance parameter to be cal- 
culated from the amount of error in a given design or measured variable. These error 
equations apply essentially to flow along streamlines or along individual blade elements. 
The errors due to more than one source can often be combined and treated as a unit. 
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For random errors, the uncertainty interval is used to indicate the combined effect of 
a number of errors on a given performance parameter. 

A number of uses of the error analysis procedures are pointed out in this report. 
The error analysis is most often used to determine the accuracy of the performance 
parameter values as an aid in interpreting design or measured data. The error analysis 
is also used to determine the critical measurement, or critical step in a design proce- 
dure, which introduces the greatest amount of error. 

The sensitivity of most of the performance parameters to design errors increased 
with increasing stagger angle or decreasing loading. Inducers, which have a combina- 
tion of high stagger angle and light loading, are especially susceptible to design errors. 
No regular pattern was noted in the trends of the sensitivity to measurement errors with 
changes in stagger angle and loading. 

The error equations, as written, apply to rotors. Modifications which allow them 
to be applied to stators were pointed out. 

A limitation of an error analysis of this type is that it depends on estimates of the 
magnitudes of design errors or measurement errors. The calculated values of errors 
in the performance parameters are only as accurate as the input error estimates. 

For pump rotor test measurements, reasonably good estimates can usually be made 
of the inherent accuracies of instrumentation such as pressure transducers and data 
recording equipment. However, errors due to probe responses to pressure fluctuations, 
flow fluctuations, and radial pressure gradients are difficult to evaluate. If a single 
probe measurement is used, circumferential variations in pressures and velocities are 
not accounted for. There is also the question of how well a stationary probe measure- 
ment will average out the blade -to-blade variations in pressures and flow angles. 

In a design procedure, it is difficult to assess the accuracy of certain inputs, 
notably loss coefficient and deviation angle. Accuracy estimates are generally obtained 
from correlations of previous test data, but the amount of suitable axial flow pump data 
is rather limited. It is also difficult to estimate the accuracy of the outlet axial veloc- 
ity, but for a different reason. In a design system, the outlet axial velocity is calcu- 
lated from the requirements of radial equilibrium and continuity. This means that the 
error in outlet velocity at a given blade element will be affected by the errors at other 
blade elements. 

The error equations were derived primarily for axial flow pumps, but are not 
limited to axial flow. All of the design error equations and measurement error equa- 
tions are correct for mixed flow impellers and centrifugal pumps if the meridional ve- 
locity is substituted for axial velocity. However, the error equations apply to a design 
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system and test procedures which are standard for axial flow pumps, but are not always 
used for mixed-flow impellers. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, January 20, 1970, 

128-31. 


20 



APPENDIX A 


D diffusion factor 

f general performance parameter 

g gravitational acceleration constant 

H total head, ft; m 

AH head rise, ft; m 

H gy net positive suction head, ft; m 

i incidence angle, deg 

k cavitation number 

r radius, ft; m 

U rotor speed, ft/sec; m/sec 

V fluid velocity, ft/sec; m/sec 

x general design or measured variable 

/3 fluid angle, angle between fluid ve- 

locity and axial direction, deg 

A uncertainty interval or finite dif- 

ference 

6 deviation angle, deg 


SYMBOLS 

r\ blade element efficiency 

cp flow coefficient 

ct solidity 

head rise coefficient 
a? total-pressure loss coefficient 

Subscripts: 
i ideal 

m,n dummy variables 
s stator 

t tip 

Z axial direction 

9 tangential direction 

1 rotor inlet 

2 rotor outlet 
Superscripts: 

’ relative to rotor 
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APPENDIX B 


ERROR ANALYSIS EQUATIONS 

The equations derived for the error analysis are identified by "des" or "meas, ” 
showing either a design error equation or a measurement error equation, respectively; 
followed by symbols showing the parameters involved. Thus des - V 2 :V Z2 identifies 
the design error equation giving the error in the performance parameter V 2 for a unit 
error in the design variable The error equations are also identified by number, 

with table I giving the equation number corresponding to each possible combination of 
performance parameter and design or measured variable. Table I also gives the figure 
numbers of the carpet plots corresponding to the equations. The table is complete for 
all of the listed performance parameters and design variables or measured variables, 
so the omission of a possible combination of variables means that the performance pa- 
rameter involved is not a function of the input design or measured variable. 


Design Equations 

The design- error equations were derived from the following equations: 


Vj = V Z1 sec (Sj 

(Bl) 

V fll = V Zl tan ' J l 

(B2) 

= v| 1 + (U 1 -V zl tan? 1 ) 2 

(B3) 

-tan"*f 1 - tan /S, ) 
v Z 1 7 

(B4) 

V 2 = V Z2 sec & 2 

(B5) 

= tan" 1 | 2 - tan ( 3 '<\ 

Vzz y 

(B6) 
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V 02 " U 2 " V Z2 ^2 

(B7) 


V 1 - Vfa + < U 2 - V z2 tan f,'/ 

(B8) 


*1 = ^< U 2 V 02 - U fV,,l) 

(B9) 


! w(Vl) 2 
1 

(BIO) 


-3 

1! 

(Bll) 


n _ 1 _ V 2 + r 2 V <92 _ r l V 01 
Vi vviirj + r 2 ) 

(B12) 


Design -Error Equations 


The resulting differential formulas are the following: 


des - V,:V zl 


(B13) 

des - V 1 :fi 1 (fig. 1) 

f— ) = JL_ tan ^ A/3 X 
yVj/ 180 1 1 

(B14) 


des - V 2 :V Z2 (fig. 2(a)) 



n 

= - cos /3 2 (tan /3g tan 0 2 - 1) 


AV 


Z2 


V, 


Z2 


(B15) 
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des - V 2 :^2 (fig- 2(b)) 


/AVc 


-2— sin cos j 3 9 sec 2 /3 9 A/3 9 
180 ^ ^ ^ ^ 


(B16) 


des - V 2 :r 2 (fig- 2(c)) 


r AV r 


= sin /? 2 cos /? 2 (tan /? 2 + tan j3 2 ) 


Ar, 


(B17) 


des - /3i:V zl (fig. 3(a)) 


A/ 3 ; = - 


^52 cos^jSUtan /3\ + tan fi-,) — 

" 1 V Z1 


(B18) 


des - I3\:(3 1 (fig. 3(b)) 


A/3j = -cos 2 /3 j sec 2 ^ A/3j 


des - P 2 :Y Z2 (In fig. 3(a) substitute 0 2 for and 0 2 for jSj.) 


(B19) 


A0 2 = - cos 2 |3 2 (tan 0 2 + tan j8 2 ) 

7T 


AV 


Z2 


Z2 


(B20) 


des - /3 2 :/3 2 (In fig. 3(b) substitute /3 2 for j3j and /3g for y3j.) 


A/3 2 = -cos 2 /3 2 sec 2 /3 2 A/3 2 


(B21) 


des - /3 2 :r 2 (In fig. 3(a) substitute /3 2 for (3^ and 0 2 for /3j; multiply result by -1. ) 

Ar„ 


A^ 2 = cos 2 /3 2 (tan 0 2 + tan 0 2 ) 

7T 


(B22) 


24 

i i nu n in i i i mu im min iii in ii i nun i hiii rim i i mi i ii in mill 



des - rp i :V zl 


where 


^zi 

*il/ V Z1 


(B23) 


des - ^/ i :/3 1 (fig. 4(a)) 


^il = 


U 1 V 01 


U f 


r l 

— tan 
r t 


= 0 for /3^ = 0 



- sec fi-, esc /3 h A/3 

180 1 1 


1 


(B24) 


(B25) 


In figure 4(a), the answer given is not a percent error in tf/^, but the error in i/a is 
given as a percentage of 4* ip corresponding to the form of equation (B25). Similarly, 
figures 4(b) to (d) give errors in i/a as percentages of i/a 2 - 

des - ^ i :V z2 (fig. 4(b)) 






= -tan @2 c °t ^2 


AV 


Z2 


(B26) 


i2/ 


Z2 


where 


des - ippP '2 ( fi S- 4 ( c )) 


Uo Vqo r o 

^ i 2 ' — = — ^2 tan ^2 

T Z r , 


Uf 


L t 


xp i2 = i/a for /3j = 0 


(B27) 
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(B28) 



• n 

180 


sec^/3 2 cot /3g A ^2 


des - i/z^r 2 (fig. 4(d)) 




- (2 + tan ySg cot 


Ar r 


(B29) 


_ 2 


des - ^:V Z1 (Multiply values from fig. 5(a) by a up y The A ip i term is obtained from 
eq. (B23). ) 


— 2 AV 7 1 

A ,ip = Aip^ + u)(p^(tan /3’^ tan /3^ - 1) 

V Z1 


(B30) 


- 2 AV Z1 

Aip = -uocp 1 for jS 1 = 0 

1 \r 1 

V Z1 


_ 2 

des - \p:/3^ (Multiply value from fig. 5(b) by c o<p-^; use eq. (B25) for Ai/^.) 


A ip = Aip^ + -^— tan sec ^l AJ 3 ! 


(B31) 


des - ^:V z2 , /3 2 , r 2 (see eqs. (B26), (B28), and (B29)) 

Aip = Aip i 


(B32) 


des - ip:v 


Aip = - 1 <p? sec^/3* 
2 1 1 


Aw 


(B33) 


des - ^ :V Zl (In fig. 2(a) substitute /3^ for £ 2 and (3^ for /3' 2 , and multiply result 
by -2; use eq. (B23) for Ai //.. Answer is error in r? as percent of 1 - 77.) 



= i_ A\p^ + 2 cos^/3j(tan j8j tan - 1) 


AV 


Z1 


Z1 


(B34) 
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des - 77:^ (In fig. 2(b) substitute p\ for 0 2 and p 1 for p' 2 , multiply result by -2; 
use eq. (B25) for A^/^.) 


_^2L j = J- At^ + — sin 0^ cos (3^ A/3^ 

1 - v ) xf^ i 90 


des - 7?:V Z2 , P# r 2 ( see eqs ' ( B26 )> ( B28 )> and ( B29 )) 



= — 
*i 


des - 77:0; 



(B35) 


(B36) 


(B37) 


In equations (B38) and (B39) use u(r 1 + r 2 )/2rj in place of a if r 2 does not equal rj. 


des - D :V zi (see fig. 6(a)) 


AD = - 


(1 - D)cos 2 /3j(tan j3j tan j3^ - 1) + tan cos ^ 


AV 


Z1 


Z1 


(B38) 


9 AV nr I 

AD = (1 - D)cos for jS 1 = 0 


V 


Z1 


des - Dt/Sj (see fig. 6(b)) 


AD = - 


180 


(1 - D)sin fi'i cos (3^ sec 2 /3j +— sec 2 j3j cos /3^ 

2cr 


Aj3, 


(B39) 


In equations (B40) to (B45) use u(r 1 + r 2 )/2r 2 in place of a if r 2 does not equal 
r l- 

des - D:V Z2 (see fig. 6(c)) 




sin A AV Z2 


2 a 


Z2 


(B40) 
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des - D:^ ( see fig- 6(d)) 


des - D:r<, 


V 2 

— AD 

v; 


/ sec j3l\ 

— (tan /3l + -) A/3p 

180 V * 2a J * 


Vi 


1 cos /3k / r, \p.\Ar 9 

- AD = (2 tan 0 2 + tan /3g ] — - 


2a 


r n + r 


2 ^ 2 / 


= M N 


r t ^i\ Ar 2 


r l + r 2 ?2/ r 2 


where (see figs. 7(a) and (b)) 


cos j3l, 

M = 

2a 


(B41) 


(B42) 


(B43) 


N = 2 tan /3 2 + t 311 @2 


(B44) 


An alternate form for /3j = 0 is 


vj cos /3l> 

— AD = 

V’ 2a 




tan ^2 


+ tan ^2 



(B45) 


Measurement Equations 

The measurement-error equations were derived from the following equations. 

Where similar equations apply to corresponding inlet and outlet parameters (such as 
and V^)* the subscript 1 or 2 is omitted. 

V z = V cos j8 (B46) 

V 0 = V sin j3 (B47) 
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(B48) 


V’ 2 = V 2 + (U - V Q ) 2 

(B49) 

P' = tan -1 ^— - tan j8| 

(B50) 

\ V Z / 


* i=-V (U 2 V 92 - U l V 9i) 

TT“ 

(B51) 


rp = 


g ah 






_ 2(if/- - wri 

U) = — - 

(vp 2 


D = 



r 2 V 02 - r l V 01 
aV’(ri + r 2 ) 


k = 


2gH sv - V 1 

(v ;> 2 


(B52) 


(B53) 


(B54) 


(B55) 


(B56) 


Measurement Error Equations 

The resulting differential formulas are 
meas-«p 1)2 :V 1>2 
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I 


I 


II 


I I ■mi ■■■■ i bi 


A cp _ AV 

cp V 


meas - <p^ 2 ( In fi £- 1 multiply result by -1. ) 



- tan P A/3 
180 


meas - cp^ 2 :U t 


,<p)~ u t 


(B57) 


(B58) 


(B59) 


meas - /3^ 2 :Y l, 2 ^ fig ’ 


A/3’ = - cos 2 /3'(tan /3' + tan /3) AY (B60) 

n V 


meas - jSj 2 :0j 2 ( In fi S- 2(a) substitute P' for 0 2 and p for /3 2 ; multiply result by 


A/3' = cos 2 /3'(tan p' tan p - 1)A p (B61) 


meas - /3j 2 :rj 2 (In fig. 3(a)) multiply by -1; subtract r(d/3'/dr) where dp'/dr is de- 
termined graphically. ) 


meas - j3j 2 :U^. 


Ar 
r 

(In fig. 3(a) multiply result by -1.) 


A/3’ = - 


d/3' 

dr 


JYY cos 2 ^'(tan P' + tan 0) 


meas - 


AIL. 


A p' = cos 2 /3’(tan p' + tan p) - 


Ux 



AU. 


U 


t 


(B62) 


(B63) 


(B64) 
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meas - ^/:AH 


a(i® {B65) 

\ 1 >) AH 


In the following equations, and are as defined in equations (B24) and (B27). 
meas - 



(B66) 


meas - 


meas - 


meas - ^/.r/Sg 


M * i = -——cp, AjSj 
1 180 r t 1 1 



= — — <P 2 Ai3 2 
1 180 r t 1 1 


(B67) 


(B68) 


(B69) 


meas - iJ/.:rj (In fig. 4(a) multiply result by -d/Sj/dr^; add (l/(p 1 )(d<p 1 /dr 1 ). ) 

/l d<?i jf d/3 A 

AiJ/. = i + sec /3 1 esc p 1 Ar 1 

1 \<Pi dr t 180 1 1 dr jy 1 

This expression is indeterminate for a /3j of 0. If p^ = 0, 


(B70) 


ir r l d/3 l 
Aif/. = -cp 

1 180 r t 1 dr x 


Ar-, 
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meas - (In fig. 4(a) substitute 0 2 for /3j, multiply result by d/3 2 /dr 2 and sub- 

tract (l/(p 2 )(d<p 2 /dr 2 ). ) 


/ 1 C ^2 7T 

A ^i = -^i21~ T~ + ~ sec 


\<P 2 


dr. 


180 


d&>\ 

CSC )Ar 0 

2 2 dr 2 / 2 


meas - 


/#i ^A 

A^/. * -f - ]Ar 9 for /3 1 = 0 

\ dr 2 r 2 / 2 1 


= -ip i 


AIL 


U + 


meas - 17: Vj, P v V 2 , /3 2 , r 2 



meas - r?: 



meas - rj:AH 


/aA _ A(AH) 

\ V ) ah 

meas - coiVj (In fig. 2(a) substitute j8j for /S 2 , ^3^ for /3 2 , multiply result by 2(1 
and add 2 times fig. 3(b). ) 


AV . 


Aoj = -2 cos 2 /3j |(1 - a>)(tan tan j3^ - 1) + sec 2 /3-|j - 


AV 


Aw 


= - 2h) cos 2 ^ - for yS^ = 0 


(B71) 

(B72) 

(B73) 

(B74) 

(B75) 

- w), 

(B76) 
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meas - w:/3^ (In fig. 8(b) substitute for /3 2 and /3j for /3 2 ; multiply result by 

-(!-«).) 


Aw = - (1 - w)cos 2 /3j(tan ^ + tan 


(B77) 


meas - w:V 2 (fig. 8(a)) 



Aw = 2 cos' 


^2 [( tan &2 


tan (3r 


1) + sec 2 0 2 J 



meas - w:/? 2 (fig. 8(b)) 


(B78) 


V 1 — TT 9 

Aw = — cos /3g(tan /3g + tan /3 2 )A/3 2 

Trt2 90 
v 2 


meas - co:r 1 (see eqs. (B76) and (B77)) 


W \f^l + /a£V*i + is. ^sv- 

A.V 1/ dr i yA/3jydri y»2 dr-^ 


(B79) 


( B80) 


Aw = 0 for constant Vj, H gv , and p x (usually p l = 0) 


meas - corrg 


Ao) 


= 2 


cosVj/V^ 

V r 2 


dr 2 dr 2/ 


(B81) 


Aw - 





dr. 



for /3j = 0 and constant H gv 
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The term di^/dr 0 denotes a change in xj/ produced by a change in only the outlet probe 
radius, that is, g/U^(dH 2 /dr 2 ) rather than g/U^(d AH/dr 2 ). However, the radial var- 
iations in are usually small compared with those in H 2 , in which case the second 
form can be used. 


meas - cotU^. (In fig. 2(c) substitute /3^ for /3 2 , /3^ for 
add 2 cosVj/**. ) 


/3 2 , multiply result by -2u>, and 


Aco = 


P , 

-2u> sin /3j cos /3j(tan ^ + tan /3j) + 2 cos^/S^ — 

vb 


AU* 


U* 


(B82) 


Aw = 


A A \jy . 

■2oJ sin^!. + 2 cos — 
1 1 2 
^ 1 


AU, 


U, 


for jSj = 0 


meas - w:AH 


Ato = - -?i- A(AH) 
V, 2 


(B83) 


r f 


In equations (B84) and (B85) use cr(r j + r 2 )/2rj in place of a if r 2 does not equal 
meas - D:Vj (fig. 6(a)) 

AD = -1(1 - D)cos“j3i(tan &\ tan - 1) +— cos (i\ tan /3 1 1 A (B84) 


(1 - D)cos 2 /3j(tan /3^ tan - 1) +— cos /3^ tan 


2 , AV 1 
AD = (1 - D)cosVi 

1 V, 


for /3j = 0 


meas - D:/3j (fig. 9(a)) 


AD = 


180 


(1 - D)cos 2 /3i(tan j3\ + tan p.) + — cos 

11 1 2a 1 


A/3, 


(B85) 
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In equations (B86) and (B87), use [<r(r ^ + r^]/2rq in place of a if r 2 does not 
equal r^. 

meas - DiVg (fig. 9(b)) 

V’ r 9 1 1 AV ? 

— -AD = cos^0k(tan ySj, tan /S 0 - 1) + — tan 0 9 cos 0k (B86) 

Vg L 1 1 1 2a 1 Z ] V 2 


meas - D:0 2 (fig- 9(c)) 


V’ 


— AD= — fees 2 


vk 


180 


cos 2 0 2 (tan @2 + tan $2) + ~~ cos ^2 A/3 2 
- 2a 


meas - D:r^ (see eqs. (B84) and (B85)) 


AD= - 


^ AD\ dV l + / AD\ d/3 l 
^V^/drj ^A^jydrj 


Ar . 


meas - D:r 2 (eqs. (B86) and (B87)) 


AD = - 


XadX^ 


( iV 2/ dr 2 l a V dr 2 


Ar, 


(B87) 


(B88) 


(B89) 


meas - D:U t (In fig. 2(c), substitute 0^ for 0 2 , 0^ for 0 2 , and multiply result by 
(1 - D). Then substitute 0 2 for 0 2 , 0 2 for 0 2 , multiply by Vg/Vj, and subtract from 
from previous result.) 


AD = 


, vk 

(1 - D)sin 0^ cos 0^(tan 0^ + tan 0^) - — =■ sin 0 2 cos 0 2 (tan 0 2 + tan 0 2 ) 

v i 


AU^ 

U + 


(B90) 
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meas - k:Vj (See fig. 2(a). Substitute j8j for /3 2 , /3j for /3 2 , multiply result by -2k, 
and add 2 times result of fig. 3(b).) 


Ak = 2jk cos 2 /3’ 1 (tan ^ tan ^ - 1) - cos 2 /3j sec 2 /3]J - 

V i 


(B91) 


2 AV 1 

Ak = -2(1 + k)cos^/3'^ 


for /3j = 0 


meas - k:/3^ (See fig. 8(b). Substitute /3'j for 0 2 and j3 l for /3 2 , and multiply result 
by k.) 


Ak = 


k— cos 2 i3 , 1 (tan + tan /3 1 )A/3 1 
90 A 


(B92) 


Ak = k-^— sin 2/31 A/3, for /3, = 0 
180 11 1 

meas - k:r^ (See fig. 2(c). Substitute j8j for /3 2 and j3j for /3 2 , multiply by -2k, 
and subtract (l/r^)(dk/drp. ) 


Ak = - 


2k sin j3j cos /3j(tan /3j + tan 0) 


+ r- 


dk 

dr. 


~|Ar, 


(B93) 


meas - k:H c „ 

bv 


Ak = 


- 2 g 

V’ 2 


AH 


sv 


(B94) 


4 
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APPENDIX C 


APPLICATION OF THE MEASUREMENT ERROR ANALYSIS 
TO SPECIFIC TEST DATA 

Two data points were selected for use in this example: the design point of the 
highly loaded axial rotor of references 3 and 4 and a typical operating point for the 84° 
helical inducer of reference 5. The rotor design point had an overall (p of 0. 451 and 
an overall \p of 0. 391, and the inducer operating point had an overall cp of 0. 0656 
and an overall ip of 0. 117. For both of these points, the flow was noncavitating. The 
84° inducer has near the highest blade angle used in current pump designs, whereas 
the rotor of references 3 and 4 has relatively low blade angles because of the high flow 
coefficient, and also has high loading. 

The values of measurement accuracy used in this example are estimates which in- 
clude the effects of more sources of error than just the instrumentation accuracies as 
given in references 4 and 5. The possible sources of error (pressure and flow fluctua- 
tions, etc. ) are discussed in the CONCLUDING REMARKS section of this report. The 
results should give an idea of what accuracy to expect in current pump test data. The 
measurement errors assumed in this study are listed in table III. 


j. 


TABLE III. - MEASURED ERROR ESTIMATES 


[Angle, ±1°; radius, 1 percent of passage 
height. ] 


Measured variable 

Highly loaded 

84° Inducer 


rotor 

(a) 


(a) 


tft 







±ft 

±m 



Head rise 

4 

1.2 

7 

2 


(180) 

(55) 

(210) 

(64) 

Inlet pressure 

3 

.9 

. 5 

. 15 


b (240) 

b (73) 

b (70) 

b (21) 

Inlet velocity head 

1.5 

.5 

. 5 

. 15 


(49) 

(15) 

(4. 2) 

(1.3) 

Outlet velocity head 

3 

.9 

2 

.6 


(120) 

(37) 

(58) 

(18) 


a Parenthetical values are maximum readings at the 


test points. 
^Absolute. 
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Using these estimates of the measurement errors, the amounts of error in the per- 
formance parameters due to the various measurement errors were calculated. The re- 
sults are shown in figures 10 to 15 for the performance parameters \p, on, 6, and 
<p g, and the measurements of p^, Vg, jSg, AH. These figures show the perform- 

ance data with solid lines outlining the limits of probable error as calculated from the 
error equations. 

To illustrate the computations, consider the errors in 6 (or Pp) due to errors in 
Vg for the highly loaded rotor. The sensitivity of 6 to errors in Vg is given in equa- 
tion meas - P'2'^2 anc * ^S ure 3(a). This figure shows that for rg = 4. 26 inches(10. 8 cm) 
at design flow, where Pp = 41. 7° and Pp = 49. 0°, a 1 percent error in Vp will give 
a -0. 51° error in 6. From the error of ±3 feet (0. 9 m) in 84 feet (26 m) of outlet ve- 
locity head (corresponding to a Vg of 73. 5 ft/sec (22.4 m/sec)), the error in Vg is 
1. 8 percent. This error in Vg results in a ±0. 9° error in 6. The corresponding 
limits of probable error are shown as solid lines in figure 14(a), which, for the radius 
rg = 4.26 inches (10. 8 cm), pass through 6 = 16. 75° + 0. 9° and 6 * 16. 75° - 0. 9°. 

Figures 10 to 14 show the uncertainty intervals in the performance parameters due 
to individual measurement errors. The differences in the relative magnitude of errors 
in comparing the rotor performance to the inducer performance may be due either to 
differences in the accuracy of the measurements or to differences in the sensitivity to 
measurement errors. In figure 13(a), the larger uncertainty intervals in <pp for the 
inducer (when compared on the basis of percentage of the nominal cpp value) are a re- 
sult of the poorer percentage accuracy assumed for the inducer Vg measurement. The 
other discernable differences in figures 10 to 14 can be explained by the effects of a 
change in stagger angle on the sensitivity to errors; these effects are listed in table II. 

For example, figures 10(a) and (c) show that the errors in \J/^ due to measurement 
errors in and pp are smaller for the inducer; this agrees with the predicted trend 
of less sensitivity at high stagger angles. The same trend for co appears in figures 
12(a) and (c). Further examples are the greater errors in <p p and 6 due to errors in 
Pp for the inducer, and the reduced errors in 6 due to errors in V^. 

Another item of interest is the relative importance of the different measurement 
errors in contributing to the total errors of the various performance parameters. For 
the highly loaded rotor, p^, Vg, and Pp contribute about equally to errors in 1//. 

(figs. 10(a) to (c)). But for the inducer, the influence of j3 ^ and Pp errors is dimin- 
ished, leaving the Vp measurement as the main source of error in The situation 

is similar for the performance parameters u> (figs. 12(a) to (d)) and 77 (not shown), 
where the p^, Pp, V 2 » and AH measurements contribute about equal errors for the 
conventional rotor performance, but only the V2 and AH measurements are impor- 
tant sources of error for the inducer performance. In fact, for the inducer, V2 was 
the main source of error for most of the performance parameters. To get better data, 
improving this particular measurement would be of prime importance. For the highly 
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loaded rotor, no one measurement could be singled out as being the greatest source of 
error. 

Since the errors in rotative speed are negligible, the errors in if/ (fig. 11) are 
caused almost entirely by errors in the measurement of head rise. Consequently, there 
is no change of sensitivity with blade angle. The percentage errors in if/ are about the 
same for the rotor and the inducer, although the different scales of the graphs make 
them appear different. 

For the test configurations used, there was no energy addition upstream of the 
blade inlet, so the errors in /3^ do not actually affect the if/^ values and have a smaller 
eff.ect on the u> values than shown in figure 12(a). However, under other conditions 
the relative size of the errors due to might be important, so the curves are included 
for comparison purposes. For all the performance parameters, the effect of radial 
probe positioning errors were small, and are not shown in the figures. Also, for 
nearly equal to zero, the effect of errors in Vj are negligible. 

The errors in the performance parameters due to the combined effects of all 
measurement errors are shown in figure 15. Equation (3) is used to calculate the total 
amount of error, as it is assumed that all of the measurement errors are random in 
nature. Thus the solid lines in the figures show the uncertainty intervals in the per- 
formance parameters. As errors in did not actually affect the values of i/a and 
had only a slight effect on the values of w, those particular errors were not included 
in the calculations. It is interesting to note that the errors in if/ and if/^ are about the 
same. Although the magnitude of the errors in deviation angle are smaller for the in- 
ducer, they are probably just as critical because of the greater sensitivity to errors in 
deviation angle at high blade angles, as noted in the design error equations. No other 
significant differences were found between the magnitude of errors in the inducer and in 
the axial rotor. 

An error analysis need not be restricted to the parameters considered herein. For 
the data of references 3 to 5, both the Venturi flow rate and the mass averaged parame- 
ters are amenable to error analysis techniques. The accuracy of the flow rate can be 
determined from estimates of the accuracy in the measurement of pressure across the 
Venturi. The errors in the mass averaged parameters can be determined using equa- 
tion (3) to get the combined effect of individual errors at several radial positions. Both 
are relatively simple applications of error analysis techniques. 
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(b) Outlet relative flow angle (eq. (B161). 
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Error in if/ as percent of oxp \ for 1° error in Pj Error in if/ as percent of <Jp\ for 1 percent error in V 
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Error in D for 1° error in 






















(VpV^ )A oj for 1 percent error in 



(a) Outlet velocity (eq. (B78)). 

Figure 8. - Effect of measurement errors on loss coefficient. 
















(a) Errors in pj. 



) .84 .88 .92 .96 1.00 .5 .6 .7 .8 .9 1.0 

Radius ratio 

(b) Errors in V 2 . 


Figure 12. - Uncertainty in total pressure-loss coefficient due to 
measurement errors. 





I I I I 

Highly loaded rotor; flow 
coefficient, 0. 451 


84 Helical inducer; flow 
coefficient, 0. 0656 a. 



(a) Errors in ideal head-rise coefficient. 


tb) Errors in head-rise coefficient. 



I I - 3 I 




(c) Errors in total -pressure loss coefficient. 




Id) Errors in rotor outlet flow coefficient. 



io i m i i j I _ 4 i I I I j | 

.80 .84 .88 .92 .96 1.00 .5 .6 .7 .8 .9 l.( 

Radius ratio 

(e) Errors In deviation angle. 

Figure 15. - Uncertainty in performance parameters due to combined 
effect of all meansurement errors. 


lASA-Langley, 1970 15 E-4704 






National Aeronautics and Space Administration 
Washington, D. C. 20546 


i 


OFFICIAL BUSINESS 


FIRST CLASS MAIL 



POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


04U 001 26 51 30S 70195 00903 

AIR FORCE WEAPONS LABORATORY /WL0L/ 
KIRTLAND AFB, NEW MEXICO 87L17 


ATT E. LOU BOWMAN, CHIEF, TECH. LIBRARY 


POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual ) Do Not Return 


r The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the tvidest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES : * Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 


published in a foreign language considered 
to merit NASA distribution in English. 


SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


ii ilium 



